Self-Consistent Theory of Normal-to-Superconducting Transition 

Leo Radzihovsky 
The James Franck Institute and Physics Department 
University of Chicago 
5640 South Ellis Avenue 
Chicago, IL 60637 



ABSTRACT 

I study the normal-to-superconducting (NS) transition within the Ginzburg-Landau 
(GL) model, taking into account the fluctuations in the m-component complex order pa- 
rameter ij^a and the vector potential A in the arbitrary dimension d, for any m. I find that 
the transition is of second-order and that the previous conclusion of the fiuctuation-driven 
first-order transition is an artifact of the breakdown of the e-expansion and the inaccuracy 
of the 1/m-expansion for physical values e = 1, m = 1. I compute the anomalous ri{d, m) 
exponent at the NS transition, and find 77(8, 1) ~ —0.38. In the m — 00 limit, r]{d, m) 
becomes exact and agrees with the 1/m-expansion. Near d = 4 the theory is also in good 
agreement with the perturbative e-expansion results for m > 183 and provides a sensible 
interpolation formula for arbitrary d and m. 

PACS numbers: 64.60.Fr, 74.20.D 
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Normal-to-nentra/-superfluid transition is one of the best understood second-order 
phase transitions with an unparaUeled agreement between theory, simulations and exper- 
iments. In contrast the problem of the normal-to-c/targe(i-superfluid i.e. the normal-to- 
superconducting (NS) transition is significantly more challenging. This problem was first 
studied twenty years ago by Halperin, Lubensky and Ma (HLM) with the GL model 
generalized to m complex component i/'q, superconducting order parameter. Using renor- 
malization group (rg) together with the first order expansion in e = 4 — d and 1/m to 
treat the gauge field and the order parameter critical fiuctuations, these authors found 
that the charge is a relevant operator that grows in the long wavelength limit. Although 
for an unphysically large number of order parameter components, m > 365.9/2 ^ 183, the 
Heisenberg fixed point (which controls the neutral superfiuid transition) was found to be 
unstable to a new critical point at a finite value of the charge, for physical superconductors 
(m = 1) no new perturbative critical point was found to terminate this charge instability. 
The authors interpreted these runaway rg fiows as a signal of a fiuctuation-driven first-order 
phase transition, providing a first example in which the fluctuations modify the order of the 
transition. Similar conclusions were also reached for the scalar electrodynamics, exactly 
in four dimensions in the context of quantum field theory. 

Although this interpretation is believed to be correct near d = 4, the conclusion of the 
first-order transition for the extreme type-II superconductors in d = 3, is most certainly 
suspect. Following the original work of Ref. Dasgupta and Halperin studied the 
problem on the lattice. Using duality arguments together with Monte Carlo techniques, 
they found that a 3d superconductor exhibits a second-order transition in the universality 
class of the (inverted) XY-model. Later Monte Carlo simulations in 3d further demon- 
strated that the nature of the transition changes from first- to second-order as one goes 
from the type-I to the extreme type-II superconductor. 

In high Tc superconductors the thermal fiuctuations are enhanced and lead to an in- 
crease of the critical region by several orders of magnitude as compared to the conventional 
superconductors. Unfortunately, even in these materials, the size of the critical region is 
still much too small to experimentally resolve the question of existence of fiuctuation- 
driven first-order NS transition, and therefore the nature of the transition appears to be 
an academic question. However, it is believed that the nematic-to-smectic-A (NA) tran- 
sition in liquid crystals is described by a model very similar to the GL gauge theory of 
the NS transition, Q and therefore the same conclusions apply to this system. In con- 
trast to superconductors, however, the NA transition is estimated to have a critical region 
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and the size of the fluctuation-driven first-order transition to be within the experimentaUy 
accessible range. Since the NA transition appears experimentaUy and predicted the- 
oretically H] to be continuous, I take this as a further indication of the breakdown of the 
perturbative e-expansion in 3d and question the conclusions of Ref. ^ . 

In this Letter, I reexamine the problem of the NS transition with analytical methods 
that do not rely on the perturbative expansion in e or 1/m. Using a non-perturbative 
method which amounts to solving approximate Dyson equations for arbitrary d and m, 
I find a nontrivial critical fixed point (e 7^ 0) that controls the NS transition. When the 
order parameter and the gauge field fiuctuations are taken into account the Heisenberg 
critical point (e = 0) controlling the neutral superfiuid transition is found to be unstable 
to this new critical point. I therefore show that in contrast to the previous conclusions 
based on the e-expansions, the 3d type-II superconductors undergo a second-order NS 
transition, consistent with the consensus described above. Besides being an independent 
prediction for the nature of the NS transition in 3d, corroborating the findings of Ref. 
, my approach has the advantage of working in arbitrary dimension and therefore sheds 
light on the question of how the 3d behavior is connected to the findings near d = 4. 

Within the GL description, the generalized superconductor is defined by the free- 
energy functional -Ff'i/'a, ^] of the m-complex-component superconducting order parameter 
t/jct and the electromagnetic vector potential A 
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where Vo ~ {T — Tc)/Tc, Qo = 2e/hc, and fXo is magnetic permeability of the normal metal. 
The choice of the Coulomb gauge V • A = 0, leads to few simplifications. 

I study the critical behavior of the NS transition within the self-consistent screening 
approximation (SCSA) that has previously been quite successfully applied to a variety 
of other problems. @ The approximation builds on the 1/m-expansion for general di- 



mensionality d. [|T0[ One writes downs the large m limit expressions for the renormalized 
interactions and propagators in terms of the bare ones and then replaces all the bare quan- 
tities by the renormalized ones thereby obtaining the large m limit of Dyson equations 
for the renormalized interactions and propagators. The advantage of this method is that 
in the limit m — » 00 it reduces to the exact 1/m result. Furthermore, while the straight 
1/m-expansion diverges for m — and therefore cannot be taken seriously quantitatively 



for the physical value of m = 1, the SCSA is perfectly well behaved in this limit and is 
therefore quantitatively more trustworthy for real superconductors. 

To simplify the analysis it is convenient to integrate out the gauge field, which can for- 
mally be done exactly since A appears at most quadratically Eq.(|l|) , giving Ftot/kBT = 
L [iVV'aP +ro|^/'aP + i^od^aP)' " ^ J. (x) (x) J, (x) + ^Tr logiD-j^)j . The last two 
terms are the long-range effective current-current interaction and the functional determi- 
nant generated from integration over A, respectively. Ji{x) = iqo{i^a^ ii^a " is 
the paramagnetic current, D~^^ = [— V'^/ {A-k^o) + 2 Qgi/'a ^Ij is the inverse of the gauge 
field propagator and can be read off from Eq.(|l]) . The treatment of type-I superconductors 
is relatively simple because the relevant temperature range lies well outside the critical re- 
gion, and therefore the order parameter fluctuations can be ignored. |jl|] By minimizing the 
effective free energy in the standard way [|10| I obtain a gauge field corrected mean-field 
theory describing the first-order NS transition previously found in Ref.|]I| . 

A more interesting and challenging regime is that of the type-II superconductors where 
the fluctuations in the order parameter field are strong and must be carefully taken into 
account. To treat this case I expand the free energy functional in powers of -i/^Q to quartic 
order, and because of the smallness of the order parameter near the NS transition I ignore 
the higher order corrections. I thereby obtain an effective field theory in terms of ipa. alone, 
with long-range self-interactions, described by an effective free energy Fef / [4>a\ 

E^llM = f ^*^^k){e+r,)^^^k) + l [ Uo{h,k2,p)r.{h-p)Mki)rp{k2+p)Mk2] 

(2) 

expressed in Fourier space with tpa{k) = J d'^xtpa{x)e~^'^^ and Jj^ = J d'^k/{27r)'^. The ef- 
fective long-range vertex of the quartic interaction is Uo{ki,k2,p) = Uo — fokiik2jPlj {p) /p'^ , 
where fo = IQirHaQo is the bare effective charge and P'[j{p) = Sij —piPj/p^- 

Using this effective free energy I write down the coupled Dyson equations for the 
renormalized tpct propagator G{k) and the renormalized quartic interactions u{p) and f{p) 

G-\k) =G-\k)+ I U{k,k-p,p)G{k-p) , (3a) 
Jp 

l + UoU^ip) l + foi^f{P) 

where U^ip) = m J^, G{p')G{p-p') and Uf{p) = -mPf-{p) / {d - I) /p^ j^, p[p'p{p')G{p - 
p') are the polarization bubbles. The diagrammatic version of these equations is displayed 
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in fig. |l|. I look for tlie long-wavelengtfi-limit solutions of the above integral equations 
for G{k), u{p) and f{p)- In general this can be done numerically with the simplification 
that near a critical point there are only two relevant length scales, and the correlation 
length ^, and therefore for example G~^{k) = k2k'^~^g{k^), with g{x) being the scaling 
function and kc is a constant that depends on the microscopies of the model. However, 
exactly at criticality, r = 0, the correlation length ^ diverges and the scaling function 
g{x —>■ oo) = 1. In this case k~^ is the only relevant length scale with correlation functions 
assuming even a simpler scaling form. In particular G~^{k) = /cj?/c^~^, integral equations 
above can be solved exactly, and 77 determined analytically, as I demonstrate below. 
Substituting the simplified scaling form for G~^{k) into Ilu{p) and Ilf{p) I find, 

UM = m/o(l-77/2, l-r7/2)A;-2V-4+2^ , Ufip) = m/o2(l-??/2, l-r7/2)fc-2V-4+2^ , 

(4) 

where I defined integrals, /o(a, 6) = j ,{p - p')-2«/p'-2& = r(a + 6 - d/2)T{d/2 - 
a)r(ci/2 - 6)/(47r)^/2/r(a)/r(6)/r(ci - a - 6), /,,(a,6) = j ^, p^ j{P - PT^y'^' = 
5,,/o2(a,6) + p,p,l22{a,h), /o2(a,6) = -T{a + b - d/2 - l)T{d/2 - a + l)T{d/2 - b + 

1) /2/(47r)^/Vr(a)/r(6)/r(rf - a - 6 + 2), 122(0, b)=r{a + b- d/2)V{d/2 - a)V{d/2 -b + 

2) /(47r)^/Vr(a)/r(6)/r(rf-a-6 + 2), and p is a unit vector. 

I first look at the Heisenberg critical point by setting fo = 0, which automatically leads 
to f{p) = from Eq. (|3£|) . Assuming that 4 — d > 2r], (this assumption will be satisfied by 
the solution for 77 for d < 4; for d > 4 the Gaussian fixed point result is recovered) in the 
long wavelength limit, p ^ 0, Huip) dominates over the 1 in the denominator of Eq.(|3£|) 
and the renormalized u{p) interaction reduces to a universal function. 

Using the renormalized version of the quartic interaction together with above equations 
and the scaling form for G{k) in Eq.(|3^) , I find that 77 is determined by m = /o(l — ?7/2, r] + 
(i/2 — 2)//o(l — 77/2, 1—77/2). For the physical superfiuids with m = 1 and d = ?> this result 
leads to 77 ^ 0.125. The implicit equation for r]{d, m) can also be expanded in e = 4 — d 
or 1/771 in which case I obtain rj^ = e^/(4m) (arbitrary m) and 77^ = 4/(37r^7?7) {d = 3), 
respectively. The large m limit of this result, 77^, by construction agrees exactly with the 
direct l/7n-expansion result. However, the (i — > 4 limit of SCSA, 77^, does not get the m 
dependence quite correctly when compared to the leading order in the e-expansion, where 
the result is 77 = e^(l + 777)/(8 + 2777)^. I note that this disagreement with e-expansion for 
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small m is expected from the fact that SCSA does not correctly account for triple vertex 



renormalization which are taken into account in the direct e-expansion for arbitrary m.[p!0 
I now apply the above calculations to the full problem of the NS transition. Allowing 
now for a nonzero charge, /<, 7^ 0, 1 arrive at the charged analogs of Eqs.(^. The expression 
for u{p) is the same, and f{p) in the asymptotic limit reduces to 

Substituting these screened interactions into U, defined by the renormalized version of Uq 
and then using Eq.(Pa[) I obtain. 



(7) 



kVk^-v = ^ f ( V\l-W2, 1-77/2) _ Io,\l-v/2A-^/2)hk,Pl^ip) 



Performing above integrals leads to the equation which determines ri[d,m) at the new 
superconducting critical point. 



/o(l- 77/2,77 + ^/2 -2) 

10(1-77/2,1-77/2) 
/o(l - 77/2, rj + d/2-l) + /o2(l - 77/2, 77 + d/2) - /22(1 - 77/2, ?7 + d/2) 

102(1-77/2,1-77/2) 
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The above implicit result for r]{d,m) reduces to 77^ = — 9e/(77i — 18) in the limit 
of ci — s> 4, for arbitrary m. In the regime where HLM find the NS critical fixed point 
(771 > 183), this e-expansion result is less than within 10% of their exact (to 0(e) ) value 
of fj^^^^ = —9e/m.^ It is important to note that although the complete SCSA result for 
77 (Eq.(||)) is well behaved as a function of m (see fig. |^), it breaks down at a critical value 
of rric = 18, when expanded in e. This suggests that the dissappearance of the critical 
point and the runaway rg flows for m < rric ~ 183 in the direct e-expansion of Ref. 
should be interpreted as the breakdown of the e-expansion rather than the fluctuation 
driven first-order transition. 

Expanding the result, Eq.(^ , for large 77i in powers of 1/771 I recover the results of 
HLM in this limit, obtaining 77^ = — 20/ n'^m ~ — 2.026/7n, for d = 3. It is important 
to note, however, that the direct 1/m-expansion leads to the value of 77^ = —2.026 (for 
d = 3,m = 1) that lies outside the physical range rj > 2 — d = —1. In contrast the SCSA 
approximation gives a sensible result of 77 = —0.38 for real superconductors, that is well 
within this physical range. For (i > 4, I recover the Gaussian fixed point, as expected since 
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the upper-critical dimension for the NS transition is d^c = 4. The SCSA then serves as a 
physical interpolation between large m and small e behavior. 

In conclusion, I find that the self-consistently modified large m- expansion for the 
generalized theory of the NS transition leads to a nontrivial critical point for 3d supercon- 
ductors. In contrast to the original HLM interpretation I predict a second-order transition 
for type-II superconductors in d = 3, in qualitative agreement with the work of Dasgupta 
and Halperin and with the related theoretical and experimental findings for the NA tran- 
sition. The results suggest a break down of the e-expansion below a critical value of m, 
while the actual NS transition remains continuous. 
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Figure Captions 

Fig. 1. Graphical representation of SCSA for renormalized propagator and interaction. 

Fig. 2. ■qSCSA _ 3^ f^ii curve) plotted as a function of m, showing improvement at 

small m compared to the direct 1/m-expansion result, i^m = — 20/(7r^m) (dashed 
curve) . The inset shows t^'^'CSA ^ function of m for various values of d. 
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